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ABSTRACT 

The surge and rotating stall post-instability behaviors of axial flow compres¬ 
sors are investigated from a bifurcation-theoretic perspective. A sequence of local 
and global bifurcations of the nonlinear system dynamics is uncovered. This in¬ 
cludes a previously unknown global bifurcation of a pair of large amplitude periodic 
solutions. Resulting from this bifurcation are a stable oscillation (“surge”) and an 
unstable oscillation (“antisurge”). The latter oscillation is found to have a decid¬ 
ing significance regarding the particular post-instability behavior experienced by 
the compressor. These results are used to reconstruct Greitzer’s (1976) findings 
regarding the manner in which post-instability behavior depends on system pa¬ 
rameters. Moreover, the results provide significant new insight deemed valuable in 
the prediction, analysis and control of stall instabilities in gas turbine jet engines. 
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I. Introduction 

A factor which severely limits the operating envelope of modern gas turbine 
jet engines is the fact that the axial flow compressor can become unstable when 
operated near its maximum achievable pressure rise. At such an operating condi¬ 
tion, a moderate disturbance can result in instability, or even loss, of the nominal 
operating point. In this circumstance, the compression system can assume one of 
basically two types of dynamic behavior. The first, known as surge , is a large am¬ 
plitude, periodic oscillation in which the compressor experiences a series of rapid 
flow reversals and recovery. The second, which is the more serious of the two, is 
known as rotating stall (or nonrecoverable stall ). This behavior is characterized 
by very inefficient engine operation at constant mass flow and pressure rise. A 
complete shutdown of the engine and subsequent re-start is usually required for 
recovery from this type of stall. 

In this paper we investigate, from a bifurcation-theoretic perspective, the 
surge and rotating stall post-instability dynamic behaviors of axial flow compres¬ 
sors. A sequence of local and global bifurcations of the nonlinear system dynamics 
is uncovered. This is used to reconstruct Greitzer’s [1] findings regarding the 
manner in which post-instability behavior depends on system parameters. More¬ 
over, the results provide significant new insight deemed valuable in the prediction, 
analysis and control of stall instabilities in gas turbine jet engines. 

Greitzer [1] introduced a lumped parameter, nonlinear model of axial flow 
compressor dynamics, simulations of which agreed with experimental results [2]. 
Among the key contributions of [1] was the introduction of a nondimensional 
parameter, B. whose value was found to be a determinant of the nature of post¬ 
instability compressor behavior (surge vs. rotating stall). Subsequently, others [7], 
[4] have sought to explain surge behavior by showing that Greitzer’s model and 
extensions thereof support a Hopf bifurcation, i.e., the emergence of a periodic 
solution from an equilibrium point as a parameter is quasistatically varied. The 
main assumption in the Hopf Bifurcation Theorem is that a pair of eigenvalues of 
the system linearization cross the imaginary axis transversely at a critical value of 
the parameter (see, e.g., [3]). 

In contrast, here we do not attribute surge in nonaxisymmetric axial flow 
compressors to this bifurcation. Rather, we present evidence for a bifurcation 
sequence which predicts the observations of Greitzer [1] for both surge and rotating 
stall. Hopf bifurcation, which has been noted by previous researchers, is accounted 
for in this work, but its role can be viewed is being of secondary significance. In this 
discussion, rotating stall and surge and their dependence on B are considered from 
a global bifurcation standpoint. Special care is taken in the choice of an analytical 
model for the steady-state compressor characteristic which closely matches that 
reported by Greitzer [1], Using this characteristic, we detect a global (so-called 
“cyclic fold”) bifurcation of a pair of large amplitude periodic solutions. One 
of these is a stable limit cycle and corresponds to the surge type of stall, while 
the second is unstable and is therefore unobserved in experiments. We refer to 
this latter unstable oscillation as the antisurge oscillation, and to the former as 
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the surge limit cycle. We find that the changing size and shape of the antisurge 
oscillation determines the critical value of B at which the post-instability behavior 
of the compression system switches from rotating stall to surge. 

The presence of the global bifurcation of the surge and antisurge oscillations 
noted above depends strongly on the assumed compressor characteristic. For ex¬ 
ample, the cubic characteristic taken in [4] does not support such a bifurcation. 
We have chosen a more complicated algebraic model of the characteristic which 
faithfully represents essential features of the experimentally obtained character¬ 
istics studied by Greitzer [1], [2]. The critical values of B reported here do not 
coincide exactly with those found in [1], but the reasons for this are apparent and 
further calculations, not reported in detail here, do result in agreement. Two main 
reasons for the apparent disagreement in critical B values are our use of a reduced 
order version of the model employed in [1] and of an approximate analytic model 
of the compressor characteristic. 

McCaughan [4] focused on an extension of Greitzer’s model (Moore and Gre¬ 
itzer [5]) which reflects nonaxisymmetric flow dynamics in addition to the axisym- 
metric flow dynamics assumed in the model of [1]. In the present work, the steady 
state compressor characteristic, being based on the experimental characteristics 
reported in [1], [2], is nonaxisymmetric. However, nonaxisymmetric flow dynamics 
are not accounted for here, although the results of [5] can be used to extend the 
present work in that direction. 

The paper is organized as follows. In Section II, we recall a nonlinear lumped 
parameter model from [1] for the dynamics of an axial flow compressor, and sum¬ 
marize the observations in [1] regarding dependence of the particular type of post¬ 
instability behavior on system parameters. In Section III, singular perturbation 
is used to reduce the order of the fourth order model of Section II, resulting in a 
second-order nonlinear model. The model used in our simulations for the steady- 
state compressor characteristic is also discussed in Section III. In Section IV, we 
present the bifurcation sequence which occurs as an important nondimensional 
parameter (B [1]) is quasistatically varied. Analytical and simulation evidence for 
this bifurcation sequence, and implications for post-instability behavior, are also 
given in Section IV. 

II. Lumped Parameter Nonlinear Model 

In this section, we recall the lumped parameter model of Greitzer [1] and 
summarize his observations regarding the circumstances of the appearance of the 
rotating stall and surge types of instabilities as a function of the nondimensional 
parameter B . 
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II-1 - The Model 


The fourth order lumped parameter model introduced by Greitzer [1] is given, 
in nondimensional variables, by: 


dmc 

dt 

dm t 
dt 

dAP 

dt 

dC 

dt 


Oh 

<1 

1 

o 

II 

(1) 

= (|)(A P-F) 

(2) 

= (|f)(”»c -m T ) 

(3) 

=d-)(c„~c) 

(4) 


Here, m c (resp. my) is the nondimensional compressor (resp. throttle) mass flow, 
A P is the nondimensional plenum pressure rise, and C is the nondimensional 
compressor pressure rise. The tilde notation used by Greitzer [1] for nondimen- 
sionalized quantities (e.g., rhc, A P) is not used here, due to the appearance here 
of only the nondimensional variables. Otherwise, the notation used above agrees 
with that of Greitzer. The steady state compressor pressure rise characteristic is 
denoted by C ss in (4), while F denotes the throttle pressure drop characteristic 
(both nondimensionalized). Three parameters appear in Eqs. (l)-(4): B, G and 
r. The first, J3, is proportional to rotor speed. Expressions for the parameters 
B, G , and r in terms of physical system parameters are given in [1], Note that 
Eqs. (l)-(3) are based on physical considerations and a simplified system repre¬ 
sentation, while Eq. (4) is phenomenological and is meant to model the inherent 
relaxation of compressor pressure rise to the steady state characteristic. Hence, 
the value of the time constant r is to be obtained experimentally. 

The compressor and throttle characteristics C ss and F. respectively, are of 
fundamental importance in the analysis of system (l)-(4). This is because these 
characteristics are the source of the nonlinearity of Eqs. (l)-(4). The steady 
state compressor characteristic C ss (for an axial flow compressor) is often mod¬ 
eled as in Figure 1, in which a hysteresis loop appears. Note that for such a plot, 
the nondimensional plenum pressure rise A P is not a single-valued function of 
the axial velocity parameter m c . Note also that the plot contains two (hystere¬ 
sis) branches which do not represent quasi-steady behavior, but rather indicate a 
“jump” (fast transient) between the lower and upper main branches of the char¬ 
acteristic. A quasi-steady characteristic may be obtained from that of Fig. 1 by 
replacing the two “jump branches” in the compressor characteristic with a single 
“transition branch” connecting the upper and lower branches of the characteristic, 
and smoothly connecting the remaining upper and lower branches of the charac¬ 
teristic. This is depicted in Figure 2 for the case in which the transition branch 
is positively sloped. With a positively sloped transition branch, the compressor 
characteristic is single-valued. Were the transition branch negatively sloped, the 
overall compressor characteristic would, at least for a range of values of mass flow 
rate, be triple-valued. Even a single-valued compressor characteristic such as that 
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depicted in Figure 2, however, can be consistent with the presence of a hysteresis 
loop. Hysteresis would occur as the throttle is slightly closed and then reopened. 

To facilitate computer simulation of the model, we shall in the following adopt 
such a single-valued compressor characteristic. As for the throttle characteristic 
F, a typical plot also appears in Figure 1. Note that this is a single-valued func¬ 
tion giving the throttle pressure drop F in terms of the throttle mass flow my. 
Moreover, the curve always has positive slope, and, for a variable area nozzle or 
valve type throttle, with ambient static pressure at the throttle discharge plane, 
we have the explicit representation 


F = 


(5) 


in terms of the flow-through compressor and throttle areas A c and Ay, respectively. 

II.2. Dynamic Behavior 

Next we proceed to give a brief summary of the influence of the value of D 
on the dynamic behavior of Eqs. (l)-(4), as discussed by Greitzer [1],[2]. All 
physical data used by Greitzer originated from an experimental three-stage axial 
compressor. For each case of interest, the model (l)-(4) was simulated and the 
results were compared with those obtained experimentally. 

To understand the manner in which initial conditions for the simulation were 
set up, we first investigate the conditions for equilibrium operation of Eqs. (l)-(4). 
At an equilibrium point, we have 


C — C ss = AP = F 


( 6 ) 


and 


m c - my (7) 

Therefore, an equilibrium point of (l)-(4) occurs at the intersection of the com¬ 
pressor and throttle characteristics, graphed relative to the independent variable 
m c (or, equivalently, my). This is the same as saying that the steady state oper¬ 
ating point is determined by the requirements that (i) the mass flow through the 
compressor and the throttle are the same, and (ii) the pressure rise through the 
compressor equals the pressure drop through the throttle. Figure 1 depicts this 
situation for a stable operating point near the peak of the steady state compressor 
characteristic. 

Each simulation and corresponding experiment reported by Greitzer [1],[2] 
begins from an initial condition obtained by perturbing the throttle line slightly 
to the left of its position in Figure 1, with the operating point initially near the 
peak of the compressor characteristic. This is depicted in Figure 3. (Figure 4 
depicts the effect of such a perturbation assuming a single-valued steady state 
compressor characteristic of the type shown in Figure 2.) Such a perturbation of 
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the throttle line may result, for instance, from slightly closing the throttle (cf. Eq. 
(5)). Depending on the values of system parameters, especially B, the system was 
then observed to either converge to the equilibrium point on the lower branch of 
the compressor characteristic (rotating stall), or to converge to a large amplitude 
periodic motion (surge). In Figures 3 and 4, the equilibrium on the lower branch 
of the compressor characteristic is denoted RSE, for rotating stall equilibrium. 

The general dependence on B of the mode of stall encountered is now re¬ 
viewed. For each of the following cases, we are concerned mainly with the trajec¬ 
tory starting from the initial condition identified above, namely that at the peak 
of the compressor characteristic, with the throttle line perturbed slightly to the 
left of the peak of the compressor characteristic (as in Figures 3 and 4). The four 
cases are referred to as Regimes 1,2,3,4. Analytical and simulation evidence will 
be given in Section IV for these behaviors. 

Regime 1: Small B. For small values of B , the trajectory settles on the rotat¬ 
ing stall equilibrium. The motion is “over-damped,” i.e., it doesn’t exhibit an 
oscillatory component. (This is depicted in Fig. 5.) 

Regime 2: Moderate B. At some point, although rotating stall is still the eventual 
limit of the trajectory, the transient is no longer over-damped, but is oscillatory. 

Regime 3: Large B. For larger values of R, the trajectory no longer seeks the 
rotating stall equilibrium point, but rather converges to a large amplitude limit 
cycle motion. This is known as surge. (This is depicted in Fig. 6.) 

Regime 4■' Very Large B. As B is increased further, the system still converges to a 
large amplitude, surge oscillation, but the nature and frequency of this limit cycle 
change drastically. Specifically, the limit cycle is now characterized by two time 
scales: a slow time scale while the trajectory remains near (“hugs”) the compressor 
characteristic, and a fast time scale for motion between the main (left and right) 
compressor branches. (A figure illustrating this motion is given in Section IV.2.) 

In the next section, the analytical model used in subsequent simulations is 
given. This involves order reduction of the model (l)-(4) as well as specification 
of the throttle and compressor characteristics. 

III. Second Order Nonlinear Model 

In this section, the surge and rotating stall instabilities and their dependence 
on B are obtained as by-products of a general nonlinear dynamic study of Eqs. 
(l)-(4). The first step is to re-draw the steady state compressor characteristic in 
a way which facilitates the analysis. 

III.l. Nonaxisymmetric Single-Valued Compressor 

Characteristic 

As noted above, the hysteresis portion of the steady state compressor char¬ 
acteristic as it is typically drawn (see Fig. 1) contains a deficiency. Namely, it 
includes two branches which represent fast, transient motions as opposed to steady 
state behavior. Instead, we choose to employ a single-valued function as in Figure 
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2. The particular characteristic is presented in Section IV.2 below, which contains 
computer simulations illustrating the paper’s main results. 

III.2. Second Order Reduced Model 

The primary variables to the compressor user are mass flow (m c ) and delivery 
pressure (A P) [1]. It is possible to approximate the fourth order compressor 
model (l)-(4) by a second order model in these variables alone. Although this 
approximation is not necessary for the results given below, it is useful in simplifying 
computations. Two basic assumptions are needed to perform the model order 
reduction: first, that the relaxation to the steady state compressor characteristic 
implied by Eq. (4) occurs in a very short time (i.e., r is small), and second, that 
the throttle inertial forces are negligible. 

The first of these assumptions enables us to neglect Eq. (4) and replace C in 
Eq. (1) by C ss ■ This is easily justified using Tikhonov’s Theorem from singular 
perturbation theory, since r may be viewed as a singular perturbation parameter. 
The second assumption implies 


A P = F(m T ), 

(8c) 

and hence that mj- is given by 


m T = F~ 1 (AP) 

(86) 

These facts imply that the primary variables m c , AP, are 
by the system 

governed approximately 

Aa = B(C„(m c ) - A P) 

(9) 

d £=( 1 B X*.-r-'W) 

(10) 


For the quadratic throttle characteristic (5), the system (9), (10) takes the form 


dm c 

dtt 

dAP 

dt 


= B(C ss (m c ) - A P) 

= (jX*. - ^(AP) 1 ' 2 ) 


( 11 ) 

( 12 ) 


Before investigating analytically the second order model (9), (10) (or (11), 
(12)), we should discuss the equilibrium points of the system. Consider either a 
hysteretic characteristic as in Fig. 1, or a single-valued characteristic as in Fig. 
2. Depictions of an “equivalent” throttle line relating A P and m c appear in Figs. 
1 and 2, along with the compressor characteristic. Equilibrium points are given 
by the points of intersection of the compressor characteristic and the throttle line. 
If the throttle line is moved to the left (closing the throttle), the rotating stall 
equilibrium point appears, as shown in Figs. 3 and 4, respectively. 



IV. Bifurcations and Stall Phenomena 

IV.1. Local analysis at rotating stall equilibrium point 

Denoting by superscript 0 an equilibrium value of a state variable, and by 
a prime differentiation with respect to the argument of a function, the Jacobian 
matrix of the second order reduced model (11), (12) is given by 

o ( C'ssi^c) - 1 

\ B~ 2 —QAA^Aq 1 B~ 2 {AP°)~ 1 / 2 

Singular perturbation results imply that as B —> 0 the eigenvalues of the Jacobian 
matrix are asymptotically given by 



B(C;,(m?) - 2il(AP“)>/ 2 ) + 0(B\ (14a) 

J* 

and 

At 

~ 2A c B(AP°y/ 2 + 

Each of the two quantities is both real and negative. (The former is negative 
because of the steepness of the throttle line, i.e., A c /At is large, and, concurrently, 
C( a (m°) is relatively small.) The negativity of the eigenvalues for small B is in 
agreement with the behavior discussed for Regime 1 in Section II.2 above, since it 
implies overdamped convergence of nearby trajectories to the equilibrium. Nothing 
can yet be said regarding the global behavior of trajectories. 

As B increases, the asymptotic formulae (14a), (14b) suggest that the eigen¬ 
values approach one another, meet, and break off in conjugate directions into the 
complex plane. This can be checked with the aid of the characteristic polynomial 
of the Jacobian matrix (13). It is more convenient to work with the following 
polynomial, which is the characteristic polynomial of 1/P times the matrix of Eq. 
(13): 


(146) 


A 2 +( 


Aj' 


'2A c B 2 (AP°y/ 2 


n> j_L_n _ sJyi}c) \ /1 r\ 

^K))H b2 (1 2 A c (AP) 1 / 2 1} 


The calculations are a bit involved. 

The conditions for a Hopf bifurcation to occur for a given value of B at a 
known equilibrium point are, from C ' ss , that 


Ax 

2A c B 2 (AP°y/ 2 


C'ss(m° c ) = 0, and 


AtC' ss {iK) 

2A c (APy/ 2 


(IGa) 


( 166 ) 
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IV.2. Bifurcation sequence and post-instability behavior 

In this section, the facts discussed above are pieced together to yield an overall 
picture of the global dynamics of the system (l)-(4) as a function of the nondimen- 
sional parameter B. The transition from over-damped to under-damped conver¬ 
gence to the rotating stall equilibrium as B is raised has already been explained. 
Also, we have seen that as B is increased further, the rotating stall equilibrium un¬ 
dergoes a Hopf bifurcation and simultaneously loses its stability. The conclusions 
we reach in this section are based on a combined use of analysis and simulation. In 
the simulations, the steady-state compressor characteristic is modeled using a com¬ 
plicated algebraic formula which need not be reproduced here. It suffices to note 
that this formula achieves a characteristic curve which is a good approximation to 
that reported by Greitzer [1], [2] for a specific three-stage compressor. 

We make some brief remarks on the simulations and the manner in which they 
are reported in the figures. The simulations were performed using the SIMNON 
software package. The figures illustrating simulation results also include the SIM¬ 
NON commands used to generate the simulations. In the SIMNON code, rh c and 
A P were denoted by mcd and pd, respectively. The notation mcdr and pdr occurs 
in the commands in figures generated by simultaneous simulations in both forward 
and reverse time. These are the variables m c and AP, respectively, computed in 
reverse time. 

Generally, for small values of 5, if the throttle line is perturbed so as to 
intersect the compressor characteristic only at the rotating stall equilibrium, then 
the latter equilibrium is the global attractor for the system. That is, every initial 
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condition converges eventually to a rotating stall operating condition. Figure 7 
depicts a simulation for the approximate model we used for B = 0.15. This figure 
illustrates the behavior just noted for small values of B. (This can be seen to 
hold under general conditions by using a Liapunov function argument applied 
to a one-dimensional reduced-order system obtained by viewing B as a singular 
perturbation parameter. A simple quadratic Liapunov function in the (scalar) 
state variable of the reduced model can be seen to suffice.) 

As B is raised past a value B\, two initially identical periodic solutions are 
born, one stable and the other unstable, and envelope the compressor characteris¬ 
tic. In the bifurcation literature, this is referred to alternatively as a cyclic fold 
bifurcation or a saddle-node bifurcation of periodic solutions. This is illustrated 
in Fig. 8, which depicts computer simulations of these periodic solutions for the 
hypothesized compressor characteristic. The unstable oscillation is found by sim¬ 
ulation in reverse time with a judicious choice of initial conditions. Note that 
the unstable periodic solution lies within the stable one. We refer to the stable 
periodic solution as the surge limit cycle (or simply surge), and to the unstable 
periodic solution as the antisurge oscillation (or simply antisurge). 

The antisurge oscillation forms the boundary of the domain of attraction 
of the rotating stall equilibrium point. 1 Note that, thus far, the peak of the 
compressor characteristic lies within the domain of attraction of the rotating stall 
equilibrium-agreeing with Greitzer’s [1] simulation results (Regimes 1 and 2 of 
Section II.2 above). Thus, for these low values of B , rotating stall must prevail 
as the steady state post-instability behavior. In addition, the unstable periodic 
solution forms the finite part of the boundary of the domain of attraction of the 
companion stable periodic solution. 

As B is increased further, the amplitude of the antisurge oscillation decreases 
steadily, until past a value B 2 , the nominal equilibrium at the peak of the com¬ 
pressor characteristic no longer lies within this periodic solution. Instead, it now 
lies within the domain of attraction of the stable periodic solution, and so is at¬ 
tracted to it. This is the initiation of surge (Regime 3 of Section II.2 above). Fig. 
9 gives results of a simulation showing, in particular, the shrinking of the unstable 
periodic solution with increasing B. (In Fig. 9, B = 0.45.) Fig. 10 shows two 
simulations for the perturbed system, both starting at the same initial condition 
(the pre-stall equilibrium), but with slightly differing values of B. In one of the 
simulations, B = 0.3757, while in the other, B — 0.3758. For the smaller value of 
B, the rotating stall equilibrium is approached (trajectory (a) of Fig. 10). For the 
larger value of B, the surge oscillation is approached (trajectory (b) of Fig. 10). 
The simulation of Fig. 11 shows conclusively that this critical value of B, i.e., the 
value at which rotating stall gives way to surge, is precisely the value for which 
the antisurge oscillation passes through the initial condition. 

As B is increased still further, the surge oscillation continues to envelope 

1 This is true for the reduced second-order dynamic model. The implication for 
the full-order model is that the antisurge oscillation lies on the boundary of the 
domain of attraction. 
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the compressor characteristic. In contrast, the unstable periodic solution shrinks, 
becoming vanishingly small, and finally collapses onto the rotating stall equilibrium 
point at some value B%. Indeed, 


£3 = £L pf 

In Fig. 12, which corresponds to B = 0.70, the unstable periodic solution is shown 
by simulation to have reached a rather small amplitude. The Hopf bifurcation 
occurring in the model under consideration is simply this merging of the now 
small-amplitude antisurge oscillation with the rotating stall equilibrium point. For 
B > B3, the rotating stall equilibrium point is no longer stable, and the unstable 
periodic solution no longer exists. The surge oscillation is now the global attractor 
of the system. Regime 3 of Section II.2 above is still in force. 

We now proceed to prove the existence of this surge oscillation for very large 
B, and in doing so determine its location and character. For large B, the surge 
limit cycle is an instance of a relaxation oscillation , very much akin to the van 
der Pol oscillation. In this regime the oscillation is characterized by two widely 
separated time scales. This was observed by Greitzer [1], and was referred to as 
“Regime 4” in Section II.2 above. The situation is depicted in Fig. 13, which is 
purely illustrative. The proof below is closely related to a standard construction 
of the van der Pol oscillation. See for instance Nayfeh [ 6 ]. 

In the following, we admit any throttle characteristic F with positive slope. 
Divide Eq. (9) by Eq. (10), to get 

drh c _ C ss (m c ) - AP 

6 dAP ~ m c — F~ 1 (AP) 1 j 

where e := B~ 2 . If e = 0, then 


AP = C ss {m c ), 

which states that solutions are relegated to the steady state compressor charac¬ 
teristic. We shall assume that e is very small, but nonzero. Referring to Figure 
13, consider a solution curve that starts at point Q. Since Q is off of the com¬ 
pressor characteristic, drhc I dAP is approximately +00 up to Q\ where it hits the 
compressor characteristic. Since dmc/dAP is approximately ±00 away from the 
compressor characteristic, the solution curve tends to stay on this characteristic. 
Also, from Eq. (10), and since point Q\ lies to the right of the throttle line in Fig. 
13, we will have that dAP/dt > 0 while the solution remains near the right portion 
of the characteristic. This means that the trajectory will move away from Q 1 up¬ 
ward, until it reaches the vicinity of Q 2 . There, the solution jumps horizontally to 
the left, since we will have dmc/dAP approximately —00. The solution then hits 
the compressor characteristic at Q 3 and again tends to remain on the characteris¬ 
tic. Since Q 3 lies to the left of the throttle line, Eq. (10) implies that dAP/dt < 0 
while the solution remains near the left portion of the characteristic. Hence the 
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solution moves downward along this branch, until it arrives at the vicinity of point 

Q, \ . There, it makes a horizontal rightward jump to Q 5. The solution exhibits a 
periodic clockwise motion thereafter, following the upper and lower branches of 
the compressor characteristic, with the exception of the horizontal jumps from Q 2 
to Q 3 and from Q 4 to Q 5 . 

V. Conclusion 

A bifurcation sequence which is useful in explaining the dependence of ob¬ 
served stall behavior of axial flow compressors on the nondimensional parameter 
B has been presented. The presence of a previously unknown global bifurcation 
of a periodic solution pair has been detected and found to be of fundamental im¬ 
portance in predicting the nature of post-instability compression system behavior. 
Specifically, an unstable antisurge oscillation has been found to be present for a 
large range of values of the parameter, and the position of the initial condition 
relative to this oscillation was found to be the decisive factor in determining post¬ 
instability behavior. 

Analysis of the limiting cases of the bifurcation sequence as well as computer 
simulation have been employed in verifying its validity. Further simulations, not 
reported here, indicate robustness of the conclusions to the compressor and throttle 
characteristics, within certain limits. Specifically, the results appear to hold for 
a large class of nonaxisymmetric compressor characteristics, but do not apply 
for axisymmetric models. They are thus expected to be useful in the analysis, 
interpretation of experimental data, and control design of axial flow compressors. 
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FIGURE 1. GENERIC COMPRESSOR AND 
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FIGURE 3. SHOWING EFFECT OF PERTURBED 
THROTTLE CHARACTERISTIC 
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FIGURE 2. SINGLE-VALUED COMPRESSOR 
AND THROTTLE CHARACTERISTICS 
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FIGURE 4. SHOWING EFFECT OF PERTRUBED 
THROTTLE CHARACTERISTIC 
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